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Abstract 

We study a system in which electrons in a two-dimensional electron g clS £1X6 
confined by a nonhomogeneous nuclear spin polarization. The system con- 
sists of a heterostructure that has non-zero nuclei spins. We show that in this 
system electrons can be confined into a dot region through a local nuclear 
spin polarization. The nuclear-spin-polarization-induced quantum dot has in- 
teresting properties indicating that electron energy levels are time-dependent 
because of the nuclear spin relaxation and diffusion processes. Electron confin- 
ing potential is a solution of diffusion equation with relaxation. Experimental 
investigations of the time-dependence of electron energy levels will result in 
more information about nuclear spin interactions in solids. 



I. INTRODUCTION 

The theoretical and experimental researches of quantum dots have attracted much at- 
tention in recent years [1]. Quantum dots are usually fabricated experimentally by applying 
lithographic and etching techniques to impose a lateral structure onto an otherwise two- 
dimensional electron system. Lateral structures introduce electrostatic potentials in the 
plane of the two-dimensional electron gas, which confines the electrons to a dot region. The 
energy levels of electrons in such quantum dots are fully quantized like in an atom. In 
such electrically confined quantum dots the confining potential can be well represented by 
a parabolic potential. 

Another method of low-dimensional structure fabrication consists of the application of 
spatially inhomogeneous magnetic fields. There has been proposed several alternative mag- 
netic structures subsequently realized experimentally. Among them: magnetic dots using 
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a scanning tunneling microscope lithographic technique [2], magnetic superlattices by the 
patterning of ferromagnetic materials integrated by semiconductors [3], type-II supercon- 
ducducting materials deposited on conventional heterostructures [4], and nonplanar two- 
dimensional electron gas (2DEG) systems grown by a molecular beam epitaxy [5]. Such 
systems were studied theoretically in a series of papers by different authors [6-14]. 

In the present paper we study a quantum dot system which is different from the quantum 
dot systems discussed above: (1) the electrons are confined through local nuclear spin polar- 
ization, (2) the confinement potential is inherently nonparabolic and time- dependent, it is a 
solution of the diffusion equation when considering relaxation, and (3) the dot contains elec- 
trons with only one spin direction. Such system was proposed for the first time in Ref. [15]. 
However, the properties of Nuclear-Spin-Polarization-Induced Quantum Dots (NSPIQD) 
have not been considered thus far and this is the motivation behind the present investigation. 
In our calculations we use some ideas from [16], where a nuclear-spin-polarization-induced 
quantum wire was proposed and investigated. 

Electron and nuclear spins interact via the contact hyperfine interaction. Once the 
nuclear spins are polarized, the charge carrier spins feel the effective hyperfine field, B^/, 
which lifts the spin degeneracy. The maximum nuclear field in GaAs can be as high as 
Bhf =5.3T in the limit that all nuclear spins are fully polarized [18]. This high level of nuclear 
spin polarization has been achieved experimentally. For example, the optical pumping of 
nuclear spins in 2DEG has demonstrated nuclear spin polarization on the order of 90%, [19]. 
A similarly high polarization has been created by quantum hall edge states (85%) [20]. The 
spin splitting due to such a hyperfine magnetic field is comparable to the Fermi energy of 
2DEG. It is important to note that the hyperfine field does not manifest itself magnetically 
due to the smallness of the nuclear magnetic moments. The electrons in the region where 
nuclear spins are polarized will preferably occupy the energetically more favorable states with 
the spins opposite to B^j. Furthermore, the nuclear polarization acts on the electrons as the 
effective confining potential. This effective confining potential can be used to create different 
nanostructures with polarized electrons in them. In this paper we consider a nuclear-spin- 
polarization-induced quantum dot (NSPIQD). 

The proposed system is depicted in Fig. 1. The nuclear spins are polarized homoge- 
neously along the z-axis perpendicular to the 2DEG in heterostructure by any other suit- 
able experimental method. For example, the optical nuclear spin polarization [22-24] or the 
transport polarization [25,21] can be used. The region where the nuclear spins are polarized 
is indicated by the cylinder in Fig. 1. The NSPIQD is created in the region of intersection 
of the 2DEG with the region of local nuclear spin polarization. The gate electrode below 
the 2DEG is used to control the number of electrons in the NSPIQD. Moreover, the sys- 
tem is subjected to an external magnetic field along the z-axis. The magnetic field plays 
an important role in the nuclear spin polarization process and, under specific conditions, 
increases the nuclear spin relaxation time. Assuming a small magnetic field, we can neglect 
it in our calculations, focusing on the effects caused by the confining hyperfine field. Our 
paper is organized as follows. In Sec. II we discuss the properties of nonhomogeneous nu- 
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clear spin polarization and calculate the evolution of initially-created hyperfine-field profile 
which is taken, for simplicity, in the Gaussian form. Time dependence of the electron states 
in NSPIQD is studied in Sec. III. The conclusions of this investigation are presented in Sec. 
IV. 



II. HYPERFINE-FIELD PROFILE 

Let us assume that the method of optical nuclear spin polarization is used [22-24] to 
create a NSPIQD. To pattern a nanostructure it is proposed to illuminate the sample locally 
by, for example, putting a mask on it. The usual optical technique allows one to create the 
light beams of the width of the order of the wavelength (~ 500nm). By using near fields 
optics the beam width can be sufficiently reduced (~ lOOnm). Hence a 1/im-size NSPIQD 
can be easily created by the modern experimental technique. 

There are two main mechanisms leading to the time dependence of the hyperfine field: 
the nuclear spin relaxation and the nuclear spin diffusion. We assume that the initial nuclear 
spin polarization is homogeneous in the z-direction. Then the hyperfine field evolution is 
described by the two-dimensional diffusion equation: 

accounting for the relaxation processes. Here D is the spin-diffusion coefficient, A is two- 
dimensional Laplace operator, and 7\ is the nuclear spin relaxation time [31,32]. The formal 
solution of Eq.(l) can be written as 

B hf = e~^ jG(r-r',t)B hf (r',t = 0)dr'. (2) 

Here G (r, t) = e ^ is the Green function of the diffusion equation and B h f (r', t — 0) is 
the initial hyperfine field profile. 

In this paper we consider NSPIQD having the cylindrical symmetry; that is, the hyperfine 
field Bhf is a function of r. In the simplest case, we can assume the initial condition to be 
of the Gaussian form: B h f (r, 0) = B exp (— The two parameters, d and B , define the 
half-width and the amplitude of the initial distribution of the hyperfine field, respectively. 
Then the solution of Eq. (1) is: 



B hf (r,t) = B e-*(l + ±-) Y 2d H^) , (3) 

where t = The value of to is the time it takes for the Bhf(0,t) to reduce by factor of 
two from t = due to the nuclear spin diffusion. The nuclear-spin relaxation time, Ti, in 
semiconductors at sufficiently low temperatures is rather long. It varies from several hours 
to a few minutes [23]. The available experimental values for the diffusion coefficient are 
D ~ 10 -13 cmV 1 for 75 As in bulk GaAs [33] and D = KT 14 cmV 1 in Alo.35Gao.e5 As [34]. 
For d = 1 and 5 fim taking D = 10~ 13 cm 2 s~ 1 we have t = 5 x 10 4 , 1.25 x 10 6 s. 



-1 - - 2 
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III. ENERGY SPECTRUM 



The microscopic description is based on the following Hamiltonian: 

H = A + \g*» B <7B hf (r, t) + U (z) (4) 

where m* is the electron effective mass, g* is the effective electron ^-factor {g* GaAs = —0.44), 
fiB is the Bohr magneton, a is the vector of Pauli matrices, B/j/ is given by Eq. (3) and 
U (z) is the 2DEG confining potential. We suppose, as is usually done for the 2DEG, that 
only the lowest sub-band, corresponding to the confinement in z-direction, is occupied and 
we can ignore the higher sub-bands. Thus, we omit z-dependence of the wave function 
in the following. The time scale introduced by a nuclear spin system is several orders of 
magnitude larger than the time scale of typical electron equilibration processes. In such a 
case the conduction electrons see a quasi-constant nuclear field. This simplifies calculation 
by avoiding the complications which would appear when solving the Schrodinger equation 
with the time dependence due to polarized nuclei. We take into account the electrons of 
only one spin direction (for which the effective potential is attractive). 

The one-electron eigenvalue problem with the attractive Gaussian potential (Eq. (3)) 
does not admit analytical solutions. Different approximate methods [26-30] were imple- 
mented to solve this problem. In the present paper, an analytical solution of the Schrodinger 
equation is found within the parabolic approximation of the hyperfine field [26]: 

B h f = a — br 2 (5) 

connected with Eq. (3) by the relations: 

B hf (0,t)=B hf (0,t) (6) 

and 

j o rB hf (0,t)dr = J o rB hf (0,t)dr. (7) 

Here Tq = a/b. Eq. (6) connects the depth of potentials, Eq. (7) provides equal nuclear-spin 

t t 

— -rp- — 

polarization for the two fields. From Eqs. (6,7) we obtain a = -Bqttx an d & = 



2d? 



RJ 



The energy spectrum for the parabolic potential (5) in units of E Q = 2 ^* d2 is given by 



g*fi B B e T i g*Li B B e 2T i 

Sn,m = ~ 1 + ± + )j IT^ ^ + H + ^ 

where n — 0, 1, ... and m — 0, ±1, ... . 

The exact solution of the Schrodinger equation with the Gaussian profile of the hyperfine 
field (Eq. (3)) was found numerically. Due to the cylindrical symmetry of the problem, the 
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wave function can be written as ip(p,(p) = -^e im< ^R(p). The equation for the radial part 
R(r) of wave function has a form 



ld_ d_ 

x dx dx 



m , B hf 



[x,t) 



x 2 'B hf (0,0) 



Rn,rn 0, 



(9) 



where x 



For d — 1 and 5 /im, 
3 hf (0,0) = 2.65 



rjd is the dimensionless coordinate and 7 = g ^^e^^ ■ 
taking m* = 0.067m e , we have £ = 0.57 x 10~ 3 , 0.023 x 10" 3 meV; for B h 
(50% nuclear spin polarization) and 5.3 T (100% nuclear spin polarization) corresponding 
energies are \g* /j, B B h f(0,0) = 3.4 x 10~ 2 and 6.8 x 10~ 2 meV. We have used the Shooting 
Method to solve Eq. (9), subjecting the solution to the following boundary conditions: 
R n ,rn (p — ► 0) = p' m ' and R n , m (p — ► oo) = 0. The results of the numerical calculations are 
presented below. 

The time- dependence of the electron energy levels in the NSPIQD is determined by the 
time-dependence of the confining hyperfine field. There are two characteristic times in the 
problem: the diffusion characteristic time to an d the relaxation characteristic time 7\. We 
can distinguish the diffusive regime, when t ~ to <S T±, the intermediate regime, t ~ to ~ T\ 
and the relaxation regime, t ~ T\ <ti to . Here t is the observation time. 

Figure 2 shows the time dependence of the electron energy levels for the Gaussian and 
parabolic potentials in the diffusion regime. We emphasize that the parabolic potential can 
be regarded as a good approximation of the Gaussian potential only for the ground state. 
The excited-state energy levels for the parabolic potential reveal large deviations from those 
for the Gaussian potential, which manifest in the degeneracy of states and in the shift of 
levels. This result is qualitatively similar to those obtained for 3D Gaussian and parabolic 
potential [26]. However, time-dependence of energy levels for both potentials show quite 
similar behavior. The number of energy levels in NSPIQD remains constant, whereas their 
depth decreases. From Eq. (8) it follows that in the diffusion regime the time-dependence 
of the energy levels in the parabolic potential is e nm {t) = £n '^f > . It can be shown that the 

energy levels in the Gaussian potential have the same t i me-dep e ndence . Substituting Eq. 
(3) into Eq. (9) and introducing the variable £ as x = £yl + t/t we obtain: 



1 d d m 2 



+ 1 + 



Rn,m 0. 



(10) 



thus 



The time-dependent factor, (l + j^j, appears in Eq. (10) only as a product with e r 
proving the statement. 

Figures 3 and 4 show the results obtained for the intermediate and relaxation regimes. On 
the contrary, the number of the energy levels in NSPIQD decreases in time in these regimes. 
This decrease occurs on the scale of T\. We can not explicitly obtain time-dependence of 
energy levels for the Gaussian potential in these regimes. The parabolic approximation of 
the hyperfine field serves as a good approximation again only for the ground energy level. 
The evolution of excited-state energy levels in the Gaussian and in the parabolic potentials 
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are different: the lifetimes of energy levels obtained in the case of the parabolic potential 
are shorter than in the case of the Gaussian potential. 

It is important to know the lifetime of the NSPIQD. We can consider electron states 
in the NSPIQD up to the moment when the confining potential depth is more than the 
temperature. Consequently, the lifetime t\ of the NSPIQD can be defined by the following 
condition: I 3 ^£^/(^0j _ where Jcb is the Boltzman constant and T is the temperature. 
Using Eq.(3), we calculate time ti for two limiting cases: T\ <C t and 7\ ^> t . In the first 
case (the relaxation regime), t\ ~ T\ In ^ 2 fcfr°^ • ^ n ^ e secon d case (the diffusion regime), 



ti ~ t ^gfa - 1 . Time-dependence of the half-width of NSPIQD is d{t) = dyjl + t/t . 
Let us estimate it at t — U. For T* = 30mK and B = 2.65T we have d(t) = d in the 
relaxation regime and d(t) = 3.6d in the diffusion regime. 



IV. CONCLUSIONS 

We have studied the electron energy levels of a NSPIQD created in the region of the in- 
tersection of a local nuclear spin polarization with a 2DEG. The properties of the NSPIQD 
are time-dependent because of the nuclear spin diffusion and relaxation. There are two 
characteristic time and three corresponding regimes: the diffusion regime, the intermediate 
regime and the relaxation regime. In the diffusion regime, the number of electron energy 
levels remains constant with time. In the relaxation and intermediate regime, the number 
of electron energy levels decreases with time. Time-dependence of the electron energy levels 
in the diffusion regime has a simple form. Since the characteristic relaxation time is pro- 
portional to the square of the NSPIQD radius at t — 0, it is possible to create NSPIQDs 
operating in different regimes using the same experimental setup. 

The numerical estimations allow us to conclude that the system under study can be 
realized experimentally. For a hyperfine field of just a few teslas, the experiment could be 
made at a temperature of the order of lOmK. The modern experimental technique allows 
one to create a region with local nuclear spin polarization of characteristic sizes ^ lOOnm, 
making the NSPIQD having a small size. The spectroscopy of the NSPIQD could be used to 
obtain some information about nuclear spin interactions in solids, for example, the nuclear 
spin relaxation time and the nuclear spin diffusion coefficient. 

It should be pointed out that a simplified model was used in this paper to describe the 
single electron states in the NSPIQD. We considered the influence of a nuclear spin-related 
hyperfine field on the electron states, whereas the electrons could also alter the nuclear spin 
dynamics. The well-known examples of such phenomena are the indirect long-range nuclear- 
spin interaction, electron-assisted mechanisms of nuclear spin relaxation and nuclear spin 
precession in an effective field created by the electrons [32]. Investigation of these effects is 
beyond the scope of this paper. Results of such investigations will be published elsewhere. 
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FIGURES 



Local nuclear spin polarization 




2DEG 



Gate 

FIG. 1. The geometry of the proposed experiment: the NSPIQD is created in the region of 
intersection of the 2DEG with the local nuclear spin polarization. 
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FIG. 2. Energy spectra of electrons in NSPIQD with initial half-width d = 1/mi and 
Bhf{r = 0, t = 0) = 2.65T as a function of time in the diffusion regime, Ti/io = 100. The 
black lines are the energy levels for parabolic potential labeled by quantum numbers (n, m) at the 
left. The other lines correspond to the energy levels for the Gaussian potential, the lines having a 
same color have the same quantum number m; the quantum number n is equal to for the lowest 
line of each color and increases by 1 for lines of the same color from bottom to top. 
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